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Abstract 

A novel procedure for the nonlinear superposition of two self-similar solutions of the heat conduc- 
tion equation with power-law nonlinearity is introduced. It is shown how the boundary conditions 
of the superposed state conflicts with self-similarity, rendering the nonlinearly superposed state to 
be a non-exact solution. It is argued that the nonlinearity couples with the presence of the scale 
so that the superposition in the linear case can give an exact solution. 
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I. INTRODUCTION 



Phenomena exhibiting self-similarity at different time and/or spatial scales are ubiquitous 



in nature 



u 



Self-similarity techniques exploit such symmetries for reducing the number of 



variables for describing the system. This also is the underlying idea in dimensional analysis 
2, 3|. Such techniques had much been explored in fluid phenomena J] and, more recently, 
in optics [5|. Self-similar solutions can only be constructed in the absence of scales having 
the dimensions of independent variables. As a result they are always endowed with extreme 
boundary conditions (BCs) representing the intermediate asymptotic behavior [l[ of a system 
away from the initial conditions and boundaries. 

A PDE may have more than one self-similar solution each corresponding to a specific BC. 
If the underlying equation is linear, its self-similar solutions can be superposed to obtain 
a solution satisfying a realistic boundary condition. Superposition of the solutions in the 
nonlinear case is not possible, but as is shown here, there is still an approximate symmetry 
to be exploited, leading to a procedure for nonlinear superposition of self-similar solutions. 



II. HEAT EQUATION AND SELF-SIMILAR SOLUTIONS 



As an example for illustrating the procedure, consider the nonlinear diffusion equation 

06 _ d 2 9 n+1 
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which is the dimensionless form of the heat conduction equation with power-law nonlinearity 



6|. This equation is also known as the porous medium equation [7J describing the flow of 
an isentropic gas through a porous medium. The equation is nonlinear for n ^ and we 
1 n the nonlinearity parameter. This equation has two well-known self-similar solutions 
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101 ]. The first solution, 
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satisfies the Dirichlet BC #(0, r) = 0, and the second solution, 



r) = r ™+2 
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satisfies the Neumann BC 9'(0, r) = where 6' = d6/d^. In both these solutions 
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The equation ([T]) can equivalently be replaced by two equations in conservative form. 
The first one is 

^ + ^ = ° (5) 

where 

is the Fick's law in dimensionless form, and the second one is 
where 

r(£,r) = $(£,T)£ + ^ +1 (8) 

It is clear from the Eqns.flS]) and ([7j) that, in the steady state, $ and V are the integration 
constants and solving equation (|5J) for 9 gives 

^(0 = (r-$0^, (9) 

which is the general solution in terms of $ and T. 

How do the time-dependent solutions given in Eqns. (j2]) and (J3]) look when written in 
terms of T(r) = — r(0,r) and $(t) = — $(0, r) corresponding to the time dependent case? 
Here the negative sign is to make the flux leaving the system from the left boundary positive. 
Using the divergence theorem for equation (JHJ), it is possible to see from the first solution 
given in equation (j2J) that 

•w-|jf«e do) 

declines as a power-law $(r) = — t~ q where 

a = l + — - — r (11) 

2(n + l) V ; 

while 

rw =sjT {Me (l2) 

vanishes (r(r) = 0). 

Similarly, it is possible to see from the second solution, given in equation (J3J), that T(r) = 
r _/3 where 



^ l ~—2 (13) 



while $(r) = 0. 



III. THE SOLUTIONS IN TERMS OF FLUXES 



We write the first solution given in equation (jSJ) in terms of $(t) as 

l/n 



(14) 



It is possible to show that this satisfies equation (pQ) even if $(r) is multiplied with a constant 
$o so that 

$(r) = -$ r- Q . (15) 
The second solution given in equation ([3]), can be written in terms of T(r) as 

l/n 



(16) 



^r(r)]^ ji - [r(r)]^- 1 } 

and it is possible to show that this satisfies equation ([1]) even if T(r) is defined as 

r(r) = r r^ (17) 

where To is a constant. Note that there is a "duality" T -h- — between the two solutions 
JH and (HHJ). 

The solutions for the linear case (n = 0), using lim n _> (i + An) l / n = e A , can be written, 
in terms of T and $, as 

0(Z,t) = -^(r)e-^ 4T (18) 

where $ = — $ r~ 3 / 2 and 

0(e,r) = r(r)e-« 2 / 4T (19) 

where T = Tor -1 / 2 , respectively. As the diffusion equation ([T]) is linear for n = 0, these two 
can be superposed to give 

e^,T)^[T(r)-mr))e~ e/4T , (20) 

and this satisfies the boundary condition #(0, r) = T(r). Note that the steady state solution 
given in Eqn.([9]) for the linear case (n = 0) becomes = T — £$. The term in the square 
brackets in Eqn.( [20l carries a form similar to the steady state solution but T and $ are 
time-dependent in the former case. 

In the linear case multiplying a solution with any constant can be absorbed into the 
constants $o or Tq. In the nonlinear case multiplying the solution with a constant does not 



give a solution, but it is possible to gauge <3>o or T to obtain the similar effect. In other 
words, multiplying the flux with a constant in the general case reduces to multiplying the 
solution with a constant. 

IV. NONLINEAR SUPERPOSITION 

The self-similar solutions given by equations ( I14jl and (fT6l) correspond to extreme bound- 
ary conditons: According to the former, there is a sink at £ = 0, such that all the heat 
carried to this boundary is totally absorbed. The latter solution describes the case in which 
there is a perfect insulator at £ = such that heat cannot leave the system from this bound- 
ary. In the steady state, the solution of which is given in equation ([9]), any ratio between 
the integration constants T and $ is possible allowing for intermediate BCs. This is not the 
case with the self-similar solutions which can only be constructed in the absence of a scale: 
If 9 and 9' are both finite at the boundary, then I ~ 9/9' is a length scale associated with 
the system. It is the aim of this paper to find the time dependent version of the general 
steady state solution given in equation (Q even though it may not be an exact solution. 

The self-similar solutions given in equations ( 1T41) and ( 1T61) appear to be the time- 
dependent versions of the steady state solution given in equation ([9]) with T = and $ = 0, 
respectively. We have written the solutions given by equations ffT4l) and (fT6l) in a "dual" 
form suggesting a "nonlinear superposition" in which we add up not the solutions themselves 
but r(r) and — £$(r) terms in separate corresponding parts of the solutions. This procedure 
gives 

9(£,t) = [r(r) -£$(r)]^r h _ A; n i_[r( r ) -£$( r )]^T-M (21) 

which reduces to equations (TH1) and ( TT6i) for r = and $o = 0, respectively. This expression 
is endowed with the BC 9(0, r) = [r(r)] 1 ^ n+1 ^ which can describe cases in which a fraction 
of the heat is absorbed at the boundary. 

The expression given in equation ( l2~lj) cannot satisfy the diffusion equation (TjQ) for a 
general value of n but is always a good approximate solution. In Figure Q] we compare 
the self-similar solutions with the numerical solution of the diffusion equation (CQ) for two 
different BCs. In both cases we have taken n = 7/3 and initiated the numerical solution from 
a Gaussian temperature distribution. The analytical solutions are time shifted as r — > r + 1 
exploiting the symmetry of the diffusion equation (TjQ) under this transformation. This allows 
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FIG. 1: Evolution of a Gaussian heat distribution on a semi infinite bar described by the diffusion 
equation ([!]), with snapshots taken at r = 0, r = 3, r = 9, r = 27 and r = 81 where r is time in 
units of diffusive time scale. Solid lines show the analytical solutions and dashed lines with data 
points (corresponding to the numerical grids) show the numerical solutions. Left panel corresponds 
to the case where we employed the Neumann boundary condition #'(0, r) = and is compared with 
the exact solution given in equation (]16p with r — > r + 1. Right panel corresponds to the case 
9(0, r) = [r(r)] 1 /( ?1+1 ^ and is compared with the approximate solution given in equation (|2ip . 

us to evaluate the analytical expression at t = 0. The left panel shows the case with 
6'(0,t) = where the exact solution given in equation f|T6l) is compared with the numerical 
solution. It is seen that it takes about ~ 10 diffusive time scales for the numerical solution 
to forget its initial configuration and to settle onto the self-similar solution. The right panel 
shows a similar comparison for the approximate analytical solution given in equation (|2ip 
for T = 0.1, $o — 1-0. We see that the approximate solution is also remarkably accurate 
after ~ 10 diffusive time scales. The analytical approximate solution is even more successful 
if n < 1. 
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V. DISCUSSION 



We have constructed an approximate solution of the nonlinear heat equation with power- 
law nonlinearity satisfying a more general BC, by using a "nonlinear superposition" of two 
self-similar solutions endowed with Dirichlet and Neumann BCs. We can not expect to find 
an exact solution for such a general BC introduces a length scale into the problem, thus, 
destroying the very condition for self-similarity. Requiring equation (j2ip be a solution of 
equation ([!]) reduces to the requirement that 



£ + „W + (£ + £ + 2W(£ + /j|r> (22) 



■) 

vanishes. The first and the third terms vanish by equations f fl5|) and (fTT|) . respectively. The 
second term, noting that 2 — a — (3 = k n , does not vanish but simplify to A; n ^T$. This 
vanishes for either r = or $ = corresponding to the well known self-similar solutions 
given in Eqns.(lH]) and ffT6j) . It also vanishes for the linear case, n = 0, as ko = 0. The 
expression given in equation ( 12T|) does not satisfy the diffusion equation (TI)) in general, but 
it is a very accurate approximate self-similar solution of it. 

The procedure can be applied to other second order nonlinear PDEs as long as they have 
two self-similar solutions and the equation can be written in terms of two PDEs in conser- 
vative form. The application of this procedure to the viscous evolution of a thin accretion 
disk in the gravitational field of a central star, a ubiquitous phenomena in astrophysics 111 ], 
will be published elsewhere later. 

Equation (|22|) vanishes when the nonlinearity parameter, n, vanishes (note that k = 0). 
As the "nonlinear superposition" technique we employed here reduces to ordinary summation 
for n = 0, this gives a way of seeing why summing up two solutions in the linear case yields 
a solution even when a scale is present: the nonlinearity term n couples with the finite 
scale in the problem and vanishing of either yields a solution. It is interesting to see that 
superposing two solutions by addition is a special case of a more general procedure employed 
in reaching equation fJITLj) from the solutions given by equations fTH|) and ffl6|) although this 
general procedure does not necessarily give an exact solution. 
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